Imaging dynamical chiral symmetry breaking: pion wave function on the light front 
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We project onto the light-front the pion's Poincare-covariant Bethe-Salpeter wave-function, ob- 
tained using two different approximations to the kernels of QCD's Dyson-Schwinger equations. At 
an hadronic scale both computed results are concave and significantly broader than the asymptotic 
distribution amplitude, (p'^^ (x) = 61(1 — x); e.g., the integral of ipTr{x)/(p'^^ (x) is 1.8 using the sim- 
plest kernel and 1.5 with the more sophisticated kernel. Independent of the kernels, the emergent 
phenomenon of dynamical chiral symmetry breaking is responsible for hardening the amplitude. 

PACS numbers: 12.38.Aw, 12.38.Lg, ll.lO.St, 14.40.Be 



The momentum-space wave- function for a nonrelativis- 
tic quantum mechanical system, ip(p,t), is a probability 
amplitude, such that \'ip{p,t)\'^ is a non-negative density 
which describes the probability that the system is de- 
scribed by momenta p at a given equal-time instant t. 
Although the replacement of certainty in classical me- 
chanics by probability in quantum mechanics was dis- 
turbing for some, the step to relativistic quantum field 
theory is still more confounding. Much of the additional 
difficulty owes to the loss of particle number conservation 
when this step is made: two systems with equal energies 
need not have the same particle content, because that 
is not conserved by Lorentz boosts, so that even inter- 
pretation via probability densities is typically lost. To 
exemplify: a charge radius cannot generally be defined 
via the overlap of two wave-functions because the initial 
and final states do not possess the same four-momentum 
and hence are not described by the same wave- function. 

Such difficulties may be circumvented by formulating 
a theory on the light-front because the eigen-functions 
of the light-front Hamiltonian are independent of the 
system's four- momentum [l], The light-front wave- 
function of an interacting quantum system therefore 
provides a connection between dynamical properties of 
the underlying relativistic quantum field theory and no- 
tions familiar from nonrclativistic quantum mechanics, 
ft can translate features that arise purely through the 
infinitely-many-body nature of relativistic quantum field 
theory into images whose interpretation is seemingly 
more straightforward. Naturally, that is only achieved 
if the light-front wave-function can be calculated. 

A phenomenon for which a quantum mechanical image 
would be desirable is dynamical chiral symmetry break- 
ing (DCSB). Strictly impossible in quantum mechanics 
with a finite number of degrees-of-freedom, this striking 
emergent feature of quantum chromodynamics (QCD), 
the strong-interaction part of the Standard Model, plays 



a critical role in forming the bulk of the visible matter 
in the Universe 0]. Expressed in numerous aspects of 
the spectrum and interactions of hadrons; e.g., the large 
splitting between parity partners d, [![ and the existence 
and location of a zero in some hadron form factors 0], 
DCSB has not yet been realised in the light-front formu- 
lation of quantum field theory. 

The impact of DCSB is expressed with particular force 
in properties of the pion. It is the pscudo-Goldstone 
boson that emerges when chiral symmetry is dynami- 
cally broken, so that its very existence as the lightest 
hadron is grounded in DCSB. As a corollary, numerous 
model-independent statements may be made about the 
pion's Bethe-Salpeter amplitude and its relationship to 
the dressed-quark propagator 0. Given that the pion's 
light-front valence-quark distribution amplitude (PDA) 
can be computed from these two quantities, their calcu- 
lation provides a means by which to expose DCSB in a 
wave-function with quantum mechanical characteristics. 

Consider, therefore, the following projection of the 
pion's Bethe-Salpeter wave function onto the light-front 



f^ip„{x) =trcD^2 / 5{n-q+ - x n ■ P) j^j ■ nx-^iq; P) , 

Jdq 

(1) 

where: f^^ is the pion's leptonic decay constant; the 
trace is over colour and spinor indices; /^^ is a Poincare- 
invariant regularization of the four-dimensional inte- 
gral, with A the ultraviolet regularization mass-scale; 
Z2(C, A) is the quark wave-function renormalisation con- 
stant, with C the renormalisation scale; n is a light-like 
four-vector, = 0; P is the pion's four-momentum, 
P^ ~ and n ■ P = —m.,^, with TOtt being the pion's 

mass; and the pion's Bethe-Salpeter wave-function 

xA<l\P) = S{q+)V^{q-P)S{q-), (2) 

with the Bethe-Salpeter amplitude, S the dressed 
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light-quark propagator, and q+ = q + rjP, q- = q — 
(1 — T])P, T] G [0,1]. Owing to Poincare covariance, no 
observable can legitimately depend on rj; i.e., the def- 
inition of the relative momentum. Using Eq. ([T]), one 
may show that the moments of the distribution; viz., 
(x™) := /g^ dx x"^ipTr{x), are given by 



^(n.Pr+i(a;™) =trcD^2 / (n • 757 • n P) ■ 

Jdq 

(3) 

The dressed-quark propagator may be expressed 

S{p) = -*7 -pavip^e) + <Js{p\e) , (4a) 

= i/[i^-pA{p\e)+B{p\e)]. (4b) 

and can be obtained from the gap equation [1, 0] : 

S-\p) = Z2{i-i-p + m}"^) 



+Zi g'^D^,y{p ~ q)—-f^,S{q)—T^{q,p), 



(5) 



where: D^i^ is the gluon propagator; T,^ the quark- 
gluon vertex; m'^™(A) the current-quark bare mass; and 
Zi{C,,K) the vertex renormalisation constant. 

The pion's amplitude may be obtained from the Bethe- 
Salpeter equation, a modern expression of which is ex- 
plained in Ref. [l3|. With r; = 1/2 in the Bcthc-Salpctcr 
equation, it is convenient to write the amplitude in the 
form 



r.(g; P) = 75 [iE^{q\ P) + i- PF^q; P) 

+ q-P-fqG^{q;P)+ a^^q^,P^H^{q- P)] 



(6) 



where the functions are even. Owing to DCSB and the 
axial- vector Ward-Takahashi identity, all scalar functions 
in Eq. ([6]) are nonzero Q. Moreover, in the chiral limit, 
which we subsequently employ exclusively, = and 



UEAq;0) = B{q^ 



(7) 



This Goldberger-Treiman-like identity, part of a near 
complete equivalence between the one-body and pseu- 
doscalar two-body problem in QCD, is a pointwise state- 
ment of Goldstone's theorem. The gap and Bethe- 
Salpetcr equations are key members of the set of Dyson- 
Schwinger equations (DSEs), which provide an efficacious 
tool for the study of hadron properties 0, Q . 

Significant features of (^tt (x) in Eq. ([T]) can be eluci- 
dated algebraically with a simple model before employ- 
ing numerical solutions for S{p), F^. To this end, with 
Am(s) = 1/[s + AP] and 7? = in Eqs. ©, ©, consider 

Sip) = [~tj-p + M]AMip^), (8) 



A r(i. + i) 



T^{q;P)^ij5— dzp{z)Al,{ql,), 

In J-1 



(10) 



where q±z = q - (l T z)P/2. Inserting Eqs. ©-(HOl) in 
Eq. , using a Fcynman paramctrisation to combine de- 
nominators, shifting the integration variable to isolate 



the integrations over Feynman parameters from that over 
the four-momentum q, and recognising that d^q-integral 
as the expression for /^r, one obtains 



r(2i/ + 2)r(m 



1) 



r(i/+ l)r(m + 2t/ + 2) 



(11) 



Suppose that = 0; i.e., the pion's Bethe-Salpeter am- 
plitude is independent of momentum and hence describes 
a point-particle, then Eq. (fTTj) yields 



r(2)r(m + 1) 



1 



(12) 



' r(l)r(m + 2) m.+ l' 

These are the moments of the distribution amplitude 

<^,(2;) = 1, (13) 

which is indeed that of a pointlike pion [ll| . 

Alternatively, consider v ^1. Then r^(A:^) ^ 1/fc^ for 
large relative momentum. This is the behaviour in QCD 
at ^ /iQ, where jig ~ 0.5 GcV is the dynamically 
generated gluon mass [I4I . u = 1 \n Eq. yields 



„ ^ r(4)r(m + 2) 
r(2)r(m + 4) 

These are the moments of 
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(TO + 3)(m + 2) 



^r(x) =6x(l-a;); 



(14) 



(15) 



viz., QCD's asymptotic PDA \\^. 

It is readily established that with Eqs. (|5))- ([TU| in 
Eq. ([3]) one obtains the "asymptotic" distribution asso- 
ciated with a {l/k'^Y vector-exchange interaction; viz.. 



Vtt{x) 



T{2v 



(16) 



Notably, the z-modulatcd dependence on g • P in Ec^. ()T[)|) 
is the critical factor in obtaining the results described 
here. To illustrate, if one uses 1/ = 1 but 2p{z) = 
S{1 — z) + 5{l + z), then point-particle moments, Eq. (fT2|) . 
are obtained even though r^(fc^) ~ 1/fc^ for fc^ IvP. 
There is a natural explanation. Namely, with such a form 
for p{z) one assigns equal probability to two distinct con- 
figurations: valence-quark with all the pion's momentum 
and valence-antiquark with none or antiquark with all 
the momentum and quark with none. In assigning equal 
weight to these two extreme configurations one has de- 
fined a bound-state with point-particle-like characteris- 
tics. It follows that deviations from the asymptotic dis- 
tribution may be expressed through Pu{z). 

We solve the gap and pion Bcthe-Salpeter equations 
numerically using the interaction in Ref. fl3 |. which 
preserves the one-loop renormalisation group behaviour 
of QCD and guarantees that the quark mass-function, 
Af(p2) = B{j)'^X'^)IA{p^X'^). is independent of the 
renormalisation point, which we choose to be C = 2 GeV. 
In completing the gap and Bethe-Salpeter kernels we em- 
ploy two different procedures and compare their results: 
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FIG. 1. Functions characterising the dressed quark propaga- 
tor. Upper panel. av{p^) ~ RL kernel: solution (open circles) 
and interpolation function (long-dashed curve); and DB ker- 
nel: solution (open squares) and interpolation function (solid 
curve). Lower panel. as{p^), with same legend. In the chiral 
limit at large p^, av{p^) ~ and crs(p^) ~ 



rainbow-ladder truncation (RL), the most widely used 
DSE computational scheme in hadron physics, detailed 
in App. A.l of Ref. [l5|: and the DCSB-improved ker- 
nels detailed in App. A. 2 of Ref. [l^ (DB), which are the 
most refined kernels currently available. Both schemes 
arc symmetry-preserving, and hence ensure Eq. ([7]), but 
the latter incorporates essentially nonperturbative effects 
associated with DCSB into the kernels, which are omitted 
in rainbow-ladder truncation and any stepwise improve- 
ment thereof • This kernel thereby exposes a key role 
played by the dressed-quark anomalous chromomagnetic 
moment in determining observable quantities (l6l | and, 
e.g., clarifies a causal connection between DCSB and the 
splitting between vector and axial- vector mesons . 

The solutions are obtained as matrices. Computation 
of the moments in Eq. ([3]) is cumbersome with such input, 
so we employ algebraic parametrisations of each array to 
serve as interpolations in evaluating the moments. For 
the quark propagator, we represent ay.s a-s meromorphic 
functions with no poles on the real p^-axis (l7j , a feature 
consistent with confinement Each scalar function in 
the pion's Bethe-Salpeter amplitude is expressed via a 
Nakanishi-like representation [l8|; i.e., through integrals 
like Eq. ([T0|) . with parameters fitted to that function's 
first two q ■ P Chebyshev moments. (Details are pre- 
sented in the Appendix.) The quality of the description 
is illustrated via the dressed-quark propagator in Fig.[TJ 



Using Eq. ([3]) it is now straightforward to compute ar- 
bitrarily many moments of the pion's PDA, {{x^"-)\m = 
1, . . . , mniax}: we typically employ mmax = 50. Since 
Gegenbauer polynomials of order a, {C"{2x — l)\n = 
0, . . . , oo}, are a complete orthonormal set on a; e [0, 1] 
with respect to the measure [x{l — x)]"- , = a — 1/2, 
they enable reconstruction of any function that vanishes 
at a; = 0,1. (N.B. (pTr{x) is even under a; O (1 — a;). 
It vanishes at the endpoints unless the interaction is 
momentum-independent.) We therefore write 



^^=(a;)=a:"-(l-a;)"- 



E 

2,4,... 



a'^Cf{2x-l) 



and minimise e,, 
A value of js 



(17) 



II 



2 ensures mean-{|(a"')'^=+V(2^'")'^' - 
l||m = 1, . . . , Tiijiax} < 1%- In using Gegenbauer-a poly- 
nomials we allow the PDA to differ from (p^^^ for any 
finite ^ and accelerate the procedure's convergence by 
optimising a. One may project our result onto a {C^^^}- 
basis, which is that used by other authors, but this incurs 
costs: requiring far more nonzero coefficients, {a^^^}, 
and introducing spurious oscillations that are typical of 
Fourier-like approximations to a simple function. 

The dashed curve in Fig. [5] is our RL result, obtained 



with Dlu = (0.87 GeV)-^ 



0.5 GeV. It is described by 



ip':'^{x) = i.74[x{i~x)r- [1+ 



Cr"(2.T-l)], (18) 



with ctRL = 0.79, = 0.0029. Projected onto a 

Gegenbauer-(a — 3/2) basis, Eq. ^TE\\ corresponds to 



,(3/2) 



0.23, 



required before a^f^^^ < 0.1 Cj'^'^'' highlights the merit 
of reconstruction via Gegenbauer-a polynomials at any 
reasonable scale, C- The merit is greater still if, as in 
lattice-QCD, one only has access to a single nontrivial 
moment. In seeking an estimate of ip^rix), it is better to 

(3 /2) 

fit a than to force a = 3/2 and infer a value for ■ 
The solid curve in Fig. [51 described by 

^°^a:) = 1.81[a:(l-a:)]"-' [1 + a^"" Cr"" {2x ^ I)] , (19) 

ttDB = 0.81, a^^ — —0.12, was obtained using the most 
sophisticated symmetry-preserving DSE kernels that are 
currently available 3, with Doj = (0.55 GeV)^, rj = 0.6. 

Projected onto a {Ci?''^}-basis, Eq. corresponds to 
a^^/^^ = 0.15. Only for j > 14 is af < 0.1 a'i^'^l 

By way of context, we note that a computation us- 
ing QCD sum rules [191 produced (p„{x = 1/2) = 1.2 ± 
0.3, which may be compared with: (p^^{l/2) = 1.16, 
ip^^{l/2) = 1.29; and with the value from the asymp- 
totic form, iy9^'^''(l/2) ~ 1.5. In addition, we find 

((2x - 1)2)RL = 0.28 , ((2a; - if)^^ = 0.25 . (20) 

Lattice-QCD ^ yields a value of 0.27 ± 0.04 for this 
moment, whereas it is 0.2 for the asymptotic distribution. 



,(3/2) 
'■14 



0.022, etc. That j > 14 is 

(3/2) 
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FIG. 2. Computed distribution amplitude at = 2 GeV. 
Curves: solid, DCSB- improved kernel (DB); dashed, rainbow- 
ladder (RL); and dotted, asymptotic distribution. 



1 1 






■ ///■' 




'•Vs 

'\V 


/// 




\\\ 


/// 







TABLE I. Representation parameters. Upper panel: Eq. (|2ip 
- the pair {x,y) represents the complex number x + iy. Lower 
panel: Eqs. (|23|l - p6|) . (Dimensioned quantities in GeV). 



RL 



Zl 



mi 



m2 



(0.44,0.014) (0.54,0.23) (0.19,0) (-1.21,-0.65) 



DB 


21 


mi 


Zs 


7712 




(0.44,0.28) 


(0.46,0.18) 


(0.12,0) (- 


1.31,-0.75) 




c' 


c" 


I/" a 


A' A" 


RL: 


E 1-cl; 


0.03 -0.71 


1.08 2.75 


1.32 1.0 




F 0.51 


c^/10 0.96 


0.0 2.78/AV 


1.09 1.0 




G 0.18 


2cp v'f 


0.0 5.73/[A'g 


f 0.94 1.0 


DB 


E 1-cl 


0.08 -0.70 


1.08 3.0 


1.41 1.0 




F 0.55 


c^/10 0.40 


0.0 3.0/Ai, 


1.13 1.0 




G -0.094 


2c'^ v'p 


0.0 1.0/[Ab 


^ 0.79 1.0 



Numerous qualitatively significant results can be read 
from Fig.[2j The most important being that DCSB is ex- 
pressed in the PDA through a marked broadening with 
respect to ^p^^ . This may be claimed because we have 
computed the PDA at a low renormalisation scale in the 
chiral limit, whereat the quark mass function owes en- 
tirely to DCSB; and, on the domain < < C^i the non- 
perturbative interactions responsible for DCSB produce 
significant structure in the dressed-quark's self-energy. 
The PDA is an integral of the pion's Bethe-Salpeter 
wave-function, whose pointwise behaviour is rigorously 
connected with that of the quark self-energy (see Eq. ([7]) 
and kindred Goldberger-Treiman relations Q). Hence, 
the structure of the pion's distribution amplitude at the 
hadronic scale is a pure expression of DCSB. As the scale 
is removed to extremely large values, phase space growth 
diminishes the impact of nonperturbative DCSB interac- 
tions, so that the PDA relaxes to its asymptotic form. 

Signficant, too, is the pointwise difference between the 
DB and RL results. It is readily understood, bearing in 
mind that low-m moments are most sensitive to ^p-„{x) 
in the neighbourhood of x = 1/2, whereas high-m mo- 
ments are sensitive to its endpoint behaviour. RL-kernels 
ignore DCSB in the quark-gluon vertex. Therefore, to 
describe a given body of phenomena, they must shift 
all DCSB-strength into the infrared behaviour of the 
quark propagator, whilst nevertheless maintaining per- 
turbative behaviour for > C,^ . This requires BijP') to 
be large ai — but drop quickly, behaviour which in- 
fluences (y9^ [x) via Eq. ([7]) . The concentration of strength 
at ~ forces large values for the small-m moments, 
which translates into a broad distribution. In contrast, 
the DB-kernel builds DCSB into the quark-gluon vertex 
and its impact is therefore shared between more elements 
of a calculation. Hence a smaller value of B{p^ = 0) is 
capable of describing the same body of phenomena; and 
this self-energy need fall less rapidly in order to reach the 
common asymptotic limit. (Using Eqs. (j?]), these remarks 
become evident in Fig.[TJ) It follows that the low-m mo- 



ments are smaller and the distribution is narrower. Both 
PDAs have the same large-x behaviour because the RL 
and DB kernels agree at ultraviolet momenta. 

Notably, one should not expect to obtain agreement 
with data for a given process by using our computed 
form of ipTr{x) in the relevant lowest-order (in coupling), 
leading-twist formula. This is illustrated well via the 
7*7 — >■ TT*^ transition form factor, G^t^^*. The dashed 
curve in Fig. [5] was obtained using a RL DSE kernel in 
that class which reproduces all uncontroversial data on 
this process uM- However, when employed in the asymp- 
totic formula [ij], the result for Q'^G^t^^* is too large by 
roughly a factor of two. Plainly, subleading contributions 
are important, at least for < 10 GeV^ and probably 
on a larger domain, as also observed elsewhere [2l|, [221 ■ 

Our PDA computations unify a diverse range of phe- 
nomena. The rainbow-ladder result, e.g., connects di- 
rectly with ab initio predictions for: tttt scattering, and 
pion electromagnetic elastic and transition form factors 
Q ; and nucleon and A properties |23| . And, although use 
of DCSB-improved kernels is just beginning, our related 
prediction for the PDA links immediately with analyses 
showing that DCSB is, e.g., responsible for both a large 
dressed-quark anomalous magnetic moment [TB] and the 
splitting between parity partners in the spectrum Q . 

The pion's PDA is the closest thing in QCD to a quan- 
tum mechanical wave function for the pion. Its hardness 
at an hadronic scale is a direct expression of DCSB. 
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Appendix. Here we describe the interpolations used in 
our evaluation of the moments in Eq. ([3]) . The dressed- 
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quark propagator is represented as [l7j 



i=i 



■ p + rrij i"f ■ p + m* 



(21) 



with Stoj ^ Vj , so that ay.s £^rc mcromorphic functions 
with no poles on the real p^-axis, a feature consistent with 
confinement 0. We find that jVn = 2 is adequate. 

With relative momentum defined via 77 = 1/2, we rep- 
resent the scalar functions in Eq. ^ {!f ^ E, F, G) by 



+ a-A%{kl) 



E^{k;P) = cl I dz p,u{z) A^.^ikl) , 
F^{k;P) - l'dzp,.{z)A-pk'Al.^{kl)., 
G-(fc; P) = c^a I'dz (z) A^^ (fc^) , 



(22) 

(23) 
(24) 
(25) 
(26) 



with Aa(s) = A^Aa(s), kl = k'^ + zk ■ P, = I - as, 
a~p = 1/A'p — Of, = 1/[Aq]'^ — aQ. H{k;P) is small, 
has little impact, and is thus neglected. 



Values of the interpolation parameters that fit our 
numerical results are presented in Tables HI Those for 
the Bethe-Salpeter amplitudes were obtained through a 
least-squares fit to the Chebyshev moments 



fn{k'') = - I dx va^5(fc; P)Un{x) , (27) 



with n = 0, 2, where Un(x) is an order-n Chebyshev poly- 
nomial of the second kind. Owing to 0(4) invariance, we 
may define x — k- P/ip, with k^ ~ 1 and P ~ {0,0, p,ip). 

The strength of the interaction detailed in Ref. [l^l 



is specified by a product: Dui 



With mo 



fixed, results for properties of ground-state vector and 
fiavour-nonsinglet pseudoscalar mesons are independent 
of the value of w G [0.4, 0.6] GeV. We use a; = 0.5 GeV. 
With the RL kernel, = 0.092 GeV is obtained 
with m^^{2GeV) = 0.87 GeV and m^^ {19 GeV) 
0.80 GeV, whilst with the DB kernel it is obtained with 
m^^{2GeV) = m^^{19GeV) = 0.55 GeV. Plainly, mul- 
tiplicative renormalisability is better preserved with the 
DB kernels. In Eq. (10) of Ref. fl, the strength of the 
dressed-quark anomalous chromomagnetic moment was 
described by a value rj = 0.65. To improve numerical sta- 
bility in the interpolations described herein, we changed 
to T] = 0.6. This increases the computed value of the ai-p 
mass-splitting by less-than 15%. 
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